CLASSIFICATION OF RADIAL SOLUTIONS TO 
LIOUVILLE SYSTEMS WITH SINGULARITIES 



CHANG-SHOU LIN AND LEI ZHANG 

Abstract. Let A = (aij) nxn be a nonnegative, symmetric, irreducible 
and invertible matrix. We prove the existence and uniqueness of radial 
solutions to the following Liouville system with singularity: 

r A Ui + E"=i««Nl^e" 3 ' (3:) =0, R 2 , i = l,...,n 

{ f n2 |a;[ ft e" <Cx) dx < oo, i = l,...,n 

where /3i,...,j3 n are constants greater than —2. If all j3iS are negative 
we prove that all solutions are radial and the linearized system is non- 
degenerate. 



1. Introduction 
In this article we consider the following singular Liouville system 
Aui + E^iOiiM^to = 0, R 2 , ie I :={!,..., n}, 



(1.1) 



J R 2 {x^ e Ui ^ dx < oo, i £ I. 



where f3\ , . . , /3 n are constants greater than 2, A — ((J'ij)nxn is a constant 
matrix that satisfies 

(HI) : A is symmetric, nonnegative, irreducible and invertible. 

A is irreducible means there is no disjoint partition of I into I\ and I2 such 
that ciij = for all i € I\ and j £ h- For the system (jl.ip . the irreducibility 
of A means (jl.ip can not be written as two independent subsystems. If 
n = 1 and an = 1, the singular Liouville system is reduced to the following 
single Liouville equation: 

(1.2) Au+ \xfe u = 0, M 2 , / \xfe u < 00. 



Prajapat-Tarantello |46j classified all the solutions to (|1.2j) and proved, on 
one hand, that if /3/2 ^ N, all solutions are radial and can be written as 

u(aQ=log , / M>°- 

V 1 8(/3/2+l) 2 l X l ^ 
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On the other hand, a solution may not be symmetric around any point if 
(3/2 € N. The proof of Prajapat-Tarantello uses properties of integrable 
system. However, the Liouville system is not integrable and we have to 
apply new methods. The purpose of this paper is to prove a classification 
theorem for all the radial solutions to (11. ip . Let u = (ui, ...,u n ) be a 
solution to (II. ip we use a = (<7i, ...,cr n ) to denote its energy: 

1 f°° 

(1.3) <n = — \x\^e u ^dx, i & I := {1, n} 

2?r Jo 

and we set A/(<t) as 

Aj(o-) = 2^2(Pi + 2)crj - ^ aijOiOy 

For Jcl, Aj(o~) is understood similarly. The main theorem of this article 
is 

Theorem 1.1. Let A satisfy (HI), /3\,...,/3 n > —2 be constants, 

(1) If u = (u\, ...,u n ) is a radial solution to then 

(1.4) A/((r) = 0, Aj(o-) > V0 C J C /. 

(2) For eoc/i a = (a\, ...,a n ) satisfying fal.4\ ), there exists a global radial 
solution u whose energy is a. 

(3) // u and v are both radial solutions to with 

poo f'OO 

Jo Jo 

Then n,(r) = Vi(5r) + (2 + /%) log 5 for some 5 > and aZZ i S 

System (ll.ip is reduced to the following form if f3\ = ... = f3 n = 0, 
f Au t + J2 je i a ij eUj =0. R2 ' 

(1.5) 

I Jm^ <oo, M 2 . 

Under the assumption (-ff 1) on ^4, a standard moving-plane argument shows 
that all Ux,..,u n are radially symmetric with respect to a common point 
(see |22j for the proof). The classification of all solutions to (|1.5p has been 
completed through the works of Chipot-Shafrir-Wolansky \22\ [23] and the 
authors [40] . Among other things Chipot et. al. prove that 

Theorem A: (Chipot-Shafrir-Wolansky) Suppose A satisfies (HI), for any 
solution u = (ui,...,u n ) to il.5}) . its energy a = (a%, a n ) belongs to the 
hypersurface 

r := {a = (o-i, ...,a n ); Aj(a) > 0, Aj(a) > 0, V0 C J C /. } 

On the other hand, for any a € T, there is a solution u of lil.5\) whose energy 
is a. 
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It can be readily verified that the energy of a solution of (|1.5p is invariant 
under rigid translations and appropriate scalings: Let n be a global solution 
to (ll.5p . then v = ■■■,v n ) defined by 

(1.6) Vi{y) = Ui (5y + xo) + 2 log 5, iel 

for any xq € R 2 and any 5 > clearly satisfies L 2 e 1 ^ = L 2 e Mi for all i £ I. 
It turns out that for any a G T, all the global solutions that have the energy 
a are related by a translation and a scaling described in fll .6jl : 

Theorem B: Suppose A satisfies (HI). Let u = (u\, ...,«„) and 

v = (v±, ...,v n ) global solutions to U.5\) such that J R2 e Ui = L2 z" 1 f or cdl 
1 G I, then v and u are related by il.6\) for some 5 > and xq gK 2 . 

Theorem A and Theorem B together give a classification of all the solu- 
tions to (jl.5p . One obvious question that Theorem 11.11 raises is, for what 
Pi, .., (3 n do all the solutions to (jl.ip have to be radially symmetric? We give 
an affirmative answer for the case of non-positive (3. 

Theorem 1.2. Let u be a solution to ( COP , A satisfy (HI). Suppose (3i € 
(—2, 0] for i € / and are not all equal to 0. Then all components of u are 
radial functions. 

Systems (jl.ip and (|1.5p and their reductions appear in many disciplines 
of mathematics and have profound background in Physics, Chemistry and 
Ecology. When fjl .5 j) is reduced to one equation, it becomes the classical 
Liouville equation 

-Aw = e u , 

which is related to finding a metric with constant Gauss curvature. In 
Physics, the Liouville equation represents the electric potential induced by 
the charge carrier in electrolytes theory [48J and is closely related to the 
abelian model in the Chern-Simons theories |29} [30l 128] . 

The Liouville systems (|1.5p(|l.ip are used to describe models in the theory 
of chemotaxis [211 E2], in the physics of charged particle beams [3 [Ml [33] . 
and in the theory of semi-conductors [22]. For applications of Liouville 
systems, see [111 [22"1 HOI I2T] and the references therein. Here we note that 
Liouville systems with singularities are of special importance in Physics and 
Geometry. For example, the single equation (jl.2p appeared in [44J as a 
limiting equation in the blow-up analysis of periodic vortices for the Chern- 
Simons theory of Jackiw and Weinberg [27] and Hong et. al. [26J. In 
geometry (jl.ip is related to finding metric with conic singularities [12\ [TS], 

ElEE]. 

It is well known that classification theorems are closely related to blowup 
analysis and degree-counting theorems. For many equations the asymptotic 
behavior of blowup solutions are approximated by global solutions. For 
example, for the Liouville equation 

An + Ve u = 0, Qcl 2 , 
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if V is a positive smooth function, blowup solutions near a blowup point can 
be well approximated by global solutions to 

Au + e u = 0, R 2 

see [El [Ml US E7] . If V is nonnegative and the blowup point happens to be 
a zero of V, the profile of blowup solutions is similar to that of the global 
solutions of (II, 2p . see [MHUISH]. We expect Theorem II .11 to be useful in the 
study of singular Liouville systems defined on Riemann surfaces or domains 
in R 2 . 

The proof of the uniqueness part of Theorem 11.11 (the third statement) is 
motivated by the authors' previous work [40J on the Liouville system with 
no singularity. The existence part (the second statement) is based on the 
uniqueness result and is therefore significantly different from the duality 
method used by Chipot. et. al. in [22j . The first statement in Theorem ll.il 
is similar to the corresponding case in |22j . 

For many applications, especially on the construction of bubbling solu- 
tions it is important to study the nondegeneracy of the linearized system. 
Our next result is concerned with the case when j5 is non-positive. 

Theorem 1.3. Let Pi £ (—2,0) for alii £ I, u = (u\,...,u n ) solve {ZHP 
corresponding to f3 = (/3i, .., f3 n ). Let (p = (0x, ■■,4>n) ^e a bounded solution 
to 

Then there exists C € M such that 4>i{r) = C(ru^(r) + 2 + /%) for all i € I. 

Remark 1.1. By Theorem \l.!% u is radial in Theorem \l.S\ 

The organization of this paper is as follows. In section two we list standard 
tools to be used in the proof of Theorem ll.il Then in section three we prove 
the three statements of Theorem 11.11 Theorem 11.21 and Theorem 11.31 are 
proved in section four and section five, respectively. Finally in the appendix 
we provide proofs for the tools used in the proof of Theorem 11.11 

Acknowledgement: Part of the work was finished when the second author 
was visiting Taida Institute for Mathematics Sciences (TIMS) in March 2011. 
He would like to thank TIMS for the warm hospitality. The second author 
is also partially supported by a grant from National Science Foundation. 

2. Preliminary results 

In this section we list a few ODE lemmas to be used in the proof of 
Theorem 11.11 Since these lemmas are standard we put their proofs in the 
appendix, in order not to disturb the main part of the paper. 

Lemma 2.1. Let u = (ui,...,u n ) be a solution to fl-l\) where A satisfies 
(HI). Then 

Ui(x) = — mi log \x\ + Cj + o(|x| _l5 ), i € /, \x\ > 1, 
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Vui(x) = -rriix/\x\ + 0(\x\ 1 ), i e I, \x\ > 1. 
w/iere 

n 

mj = ^ aijfTj > 2 + Pi, i e I, 

3=1 

poo n 

Ci = Ui (0)+ / log rV aij r fr +1 e u ^dr 
5 is some positive small number. 

Remark 2.1. u is not assumed to be radial in Lemma \2.1\ 

The next lemma is on the linearized system of (jl.ip expanded along a 
radial solution u: 

(2.1) (r$)' + Yj "•J r ''' ( "' r °j = °> * e J - 

i 

Lemma 2.2. Let 4> = (0i> ■■■,4>n) satisfy \2.1\) with (3$ > —2 /or a// i £ I, 

then <fri(r) = O(logr) at infinity for i £ /. 

Lemma 2.3. Lei A satisfy (HI), /3j > —2 for i € /, t/ien /or any c\, c n € 

R, there is a unique solution to 



I'M + ^-(r) + E"=i ayr&e^W =0, i = 1, ..,n, 
Ui(0) = Ci, i = 1, ..,n 



(2.2) 

i/iai exists for all r > 0. 

Remark 2.2. u may not have finite energy. 

If we further know that an > for all i, then the solution has a finite 
energy: 

Lemma 2.4. Let an > for all i € L, then for all c\, ...,c n € M., there exists 
a solution to 

' < + M( r ) + Ejej Oi^e^M = 0, < r < oo, » e I, 

(2.3) < / °° e^W^+'dr < oo, i G I, 

k Ui(0) = Q, i£ I. 
Lemma 2.5. Let (p be a solution of 

(r<Pi(r)Y + Y%=i Cj'- 1 - ' l <" ■(:>,(!■) = 0, < r < oo, 

&(0) = 0, Vie I. 
Suppose Pi > —2 for all i £ I, then 4>i = /or aZZ i £ L. 
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3. Proof of Theorem ll.lt 

3.1. The proof of the first statement of Theorem 11.11 

Lemma 3.1. Let u = (ui,..,u n ) be a radial solution of M.l\) with A satis- 
fying (HI) and /3, > —2 for all i. Then 

A/(<r) = 0, Aj(a) > 0, V0 C J C /. 

Proof of Lemma I3.lt This proof uses the same idea as in [22J. Let 

Ui(t) = Ui(e*), then 

(3.1) ^(t) — > —mi as t — > oo. 
The equation for Uj(i) is 

n 

(3.2) «J'(t) + ^ Oij-e^^W =0, t g R, iel. 

3=1 

Let Zi(t) = Y^j=i CL tJ Uj(t), then z[(t) — > —Oi as t — > oo. (|3.2p can be rewritten 
as 

(3.3) z'lit) = _ e ( 2 +ft)*+Ei=i^^, i G /. 

Clearly (i) < for all i G I and all t G R. Let u>j(i) = z-(t), then by (pOj) 
and (ptTD 

oo) = 0, tUj(i) < Vt, u)j(oo) = —Oi. 
In addition we have w\{— oo) = w4(oo) = 0. Using the definition of wi we 
differentiate (I3.3P to obtain 

(3.4) w'l{t) = (2 + ft)«4(t) + 

Taking the summation for i G / in (13. 4h we can write (I3.4D as 

^t«?(t) - 53(2 + A)^(t) = 53 l^M^t))'. 

iei iel i,jei 

Integrating t from — oo to oo we obtain Aj(a) = 0. For J C I, summation 
for i G J in (13.41) leads to 



£X'(f) " J> + - ^ J2 a tj (wiw 3 y(t) = Yl <HA(t)vj(t). 

ieJ ieJ i,jeJ ieJ,jei\J 

Integrating the above for t G (— oo, oo) we obtain 

1 . r°° 
^( 2 ~ 2 Y a ^°i°i = Y ai i I w'i(t)wj(t)dt. 

ieJ i,jeJ ieJ,jeJ\J 

Since the irreducibility of A means that there exists a™ > for i G J,j G J\J 
we see the right hand side of the above is strictly positive because w'^t) < 
for all i G / and t G R and wt(t) < for all i G I and i G K. Thus we have 
obtained Aj(cr) > 0. Lemma |3. II is established. □ 
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3.2. The proof of the third statement of Theorem 11.11 
Proposition 3.1. Let u and v both be radial solutions to such that 

pOO POO 

/ r ft+V i(r) dr = / r ft+ V l(r) dr, i G I. 
Jo Jo 

Then Uj(r) = Vi(5r) + (2 + /3j) log 5 for some 5 > and all i G /. 

To prove Proposition 13.11 we first establish a uniqueness result for the 
linearized system: 

Lemma 3.2. Let cp = (<pi, cf> n ) be a bounded solution of \2. then 
fa(r) = C{ru' i (r) + 2 + p i ) for all iel. 

Proof of Lemma 13.2b Let 

<jP = (ru[{r) + 2 + ..,n4(r) + 2 + /?„). 

Then by computation i^ is a solution to the linearized system. Suppose 
there exists another bounded solution cp l which is not a multiple of </>°. 
Without loss of generality we assume (fr\(0) = 0, as by Lemma [231 one of 
4>j(0) must be different from 2 + /3j. To derive a contradiction we set 

S = {a; 3 a bounded solution <f> = (cfti, n ) such that <^>i(0) = 2 + ft, 

0i(O) = < 3 + /%;, i = 2,...,re, a = min{2 + ft, a 2 , a n }, 

rr 

/ e Ul(s) (/>i(s)s 1+A ^ > 0, Vr>0, i = l,...,n. }. 

First we see that 2 + min{ft, /3 n } € S 1 . Indeed the expression of (fP gives 

[ r s 1+ ^e u * {s) (j) i{s)ds = r 2+/ V* (r) > 0. 
Jo 

Next we observe that S is a bounded set. Indeed, suppose a < is in S, 
let <j> be the function corresponding to a, then 3j G / such that </>j(0) = a. 
This leads to JT s 1+ ^ e u ^ s '(f>j(s)ds < for r small, a contradiction to the 
definition of S*. Let a be the infimum of S and let a fc = (a k , ...,«„) G S be 
a sequence in 5 that tends to a from above. Suppose cp k = (<j)\, ■ ■■,4> k l ) is 
the solution corresponding to a k , then we claim that cp k converges to cf) = 
((/>!,..., n ), which is also a bounded solution with the strict monotonicity 
property described in S. Indeed, let ifj 171 = (ip™, be the solution to 
the linearized system such that ipf 1 ^) = 5™. Then by Lemma 12.51 

n 

\ „k „i,m 
9=2^ a mW ■ 

m=l 

Here we recall that by Lemma 12.21 ibl n (r) = O(logr) for r large. Since 
a < a k < 3 + Pi for z 6 7 and all k. Along a subsequence cp k tends to (/> over 
all compact subsets of R. The monotonicity property of ct> k implies 

l e u ^ s) ^i{s)s l+ ^ds > 0, Vi G /, Vr > 0. 
Jo 
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On the other hand, since 4> k are all bounded functions, for each <p k we 
find r\ — > oo as I — > oo such that ri((f) k y(ri) ~~ > 0- From the equation for rfa? 
we have 

/ Y\ a ijr ^ +1 e u ^ r Wj(r)dr = 0, V« G /. 
•A) i=1 

Since A is invertible 

= / e ni(r) ^(r)r ft+1 ( ir = V a* / e^-i/f (r)r ft+1 dr. 

m=l ^ 

/>oo 

Since ip™-(r) = O(logr), / e Ul ( r ^ m (r)r ft+1 cir is well defined, we let a k ->■ 

Jo 

(qi, a n ) to obtain 

/■oo 

(3.5) / e u ^fa(s)s^ +1 ds = 0, Vi € /. 

Jo 

As a consequence of (I3.5j) . is bounded. Indeed, the equation for is 
(r#(r))' = -^ayr^ +1 e^W^(r), r > 0. 

Using 0j(r) = O(logr), r ft+ 2 e «iW _ for some 6 > (Lemma I2.ip 

and (|3.5p we know 

/■oo 

z u ^ s) fa{s)s^ +1 ds = - / e Ui(fl) &(s)s ft+1 ds = O^" 5 / 2 ) 



for r large. Thus fa^{r) = 0(r _1_<5 ) for all r large, which implies that fa is 
bounded. Since each fa is a non-increasing function, ([3.5|) implies that fa 
decreases to a negative constant when r — > oo. Indeed, by (|3.5p either ^ = 
or fa decreases to a negative constant. The first possibility does not exist, 
because the fact fa(0) = 2+ fa > implies that fa decreases into a negative 
constant at infinity. Also JJ" s 1+ ^ 1 e Ul ^fa(s)ds > for all r. Consequently 
for all i in the set I\ := {i G /; an > }, 

r#(r) < -oa / s 1+/3l e UlW 0i(s)ds < 0, Vr > 0. 



Therefore ^ strictly decreases to a negative constant for all i G I\. We can 
further define 

h '■= {i G ^; o«j > for some j G Ji. }. 

By the same reason as above fa decreases to a negative constant at infinity 
for all 16/2- By the irreducibility of A all the components of </> decrease to 
negative constants at infinity. 

Now we claim that a — e G S for e > small. To see this, consider cp + tcj) 1 
for \t\ sufficiently small. Recall that <f>\(0) = 0, thus fa(0)+t(p\(0) =2 + fa. 
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Clearly 4> + t<j) 1 solves (|2.ip . By choosing t positive or negative with \t\ small 
we can make 

min^(O) + t(j)}(0) = a - e > 0. 
■iei 

Since (/> + tcp 1 is bounded we have 

e u *(4>i + t(j)})s^ +1 ds = 0, i = l,...,n. 



Since (j>i{r) tends to a negative constant as r — > oo and (j) 1 is bounded, we 
know for r large and \t\ small 



e u *(<pi(s) + t<t> l i {s))s Pi+1 ds < 0. 



Consequently 



2 w<(s) (<&(s) + t(pj(s))s^ +1 ds > Vr > 0. 



Thus a — e € 5 for some e > small, a contradiction to the definition of a. 
Lemma 13.21 is established. □ 



Proof of Proposition I3.lt We shall consider 

' <0) + + Ejei aijrtoeWM = 0, < r < oo, 

/ °° e «iW r ft+i^r < oo, Vi E I, 
k iti(0) = Ci, i = 1, u n (0) = 0. 



(3.6) 



Let 

n 2 



{a = (a 1 ,..,a n ); Aj(<r) = 0, Aj(<t) > 0, V0 C J C /. }. 

LTi := {C = (ci, ..,c n _i); (|3.6p has a solution. }. 

Note that by Lemma [2T41 11 1 = M n_1 if > for all i. We claim that the 
mapping from IIi to II2 is locally one to one. Indeed, let M be the following 
matrix: 

/ d Cl at ... d Cn _ x o x 

M= 

\ 9 Cl fJ n _i ... 9 Cn _ 1 cr„_i 

We claim that M is nonsingular. We prove this claim by contradiction. 
Suppose there exists a non-zero vector D = (di, d n _i) T such that MD = 
0. Then by setting 7 = d\c\ + ... + d n -\c n -\ we have 

(3.7) <9 7 <7i = (9 7 <72 = ... = (9 7 <7 n _i = 0. 

For n 2 , A/(cr) = reads 
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By differentiating both sides with respect to 7 we have 

* 3 

Since Aj(<t) > implies ~Yl,j a ij a i > 2 + /??,, (|3.7p implies <9 7 <7 n = 0. Set 
4>i = d-yUi (i E I), then = (0i, ...,4> n ) satisfies (|2.ip and 

0i(O)=di, i = l,...,n-l, n (O)=O. 

From <9 7 <7j = (i € /) we have 

/>oo 

(3.8) / e u *<j)i(s)s 1+ ^ds = 0, i€l. 
Jo 

As a consequence of (j3.8j) . <j) is bounded. Indeed, integrating (|2. 1 j) from to 

r 

r$(r) = - [ r y^a ijS 1+ he u > (s Uj(.s)ds 
Jo j 

aij8 1+p ie^*Uj(s)d8 = 0{r- & ) 

for some (5 > 0. Therefore </>'(r) = 0(r~ 1-5 ), which proves that <fii is 
bounded. By Lemma I3T21 <fo = c(ru^ + 2 + then we see immediately that 
c = because </> n (0) = 0, this is not possible because not all dj's are zero. 
Therefore we have proved that M is nonsingular for all C = (ci, ...,c n -i) € 
IIi. 

We further assert that there is one-to-one correspondence between LTi and 
II2. This is proved in two steps as follows. 
Case 1: an > 0, i G I. 

In this case, by Lemma 12.41 IT1 = M™ -1 . The mapping from LTi to LT2 is 
proper and locally one to one. Here we claim that LT2 is simply connected. 
Assuming this, since R n_1 and LT2 are simply connected, there is one to one 
correspondence between them. Let u = (u%,...,u n ) and v = (vi,..,v n ) be 
two radial solutions such that u n (0) = v n (0) = 0, f R2 \x\^e Ul = L 2 (x^e" 1 
(i G I). Then Ui(0) = Vi(0) for i = l,...,n — 1. By Lemma 12.31 m = Vi 
for all i G I. Now we prove that LT2 is simply connected. Indeed, using 
mi = Y2j a ij a j-, A/((t) = can be written as 

(3.9) Y, a%J ( 2 + ^)( 2 + ft) = E ° 4j ( mi " 2 " &X m i " 2 - Pj). 

Therefore LT2 is part of a quadratic surface, the boundary of which is re- 
stricted by A j t [a) = where J« is / with the index i removed. A j i {a) > 
reads 

mi - 2- Pi> -^o-i- 

In another word in the coordinate system represented by rrij, we use n coor- 
dinate planes to bound the quadratic hypersurface described in (|3.9p . Other 
restrictions Aj > 0, when J is obtained from / with at least two indices 
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removed, do not affect the topological information of Aj(a) = 0. Thus II2 is 
a part of the quadratic hyper-surface in the first quadrant and is therefore 
simply connected. Proposition 13.11 is proved in this case. 

Case 2: There exists io such that ai 0t i Q = 0. We prove this case by 
a contradiction. Suppose c k = (cf, c^_ 1 ) (k = 1,2) are two distinct 
points on IIi that correspond to the same energy: let u , u 2 be two solutions 
corresponding to c 1 and c 2 respectively such that 



e u\{r) r l+pi dr= / e uf (r) r l+Pi dr = a . 

Jo 



i G /. 



(da 



1) is nonsingular at c 1 and c 2 , there 



Since the matrix J e iV1 (n-i)x(n-i) 

is a one-to-one mapping between a neighborhood of c k to a neighborhood 
of a in H2. Since c 1 7^ c 2 , we choose the neighborhoods around them to be 
disjoint. 

Now consider a perturbation system 

f <(r) + ±<(r) + E je /(«ij + e&ii)rhe"* =0, r > 0, t € I, 



(3.10) 



J °° r^ +1 e Ul dr < 00, i £ I, 
_ Ul (0) = ci, ...«„_! (0) = c^-i, u n (0) = 0. 



Let it ' e be the solution to ()3.10|) that corresponds to the initial condition 



J; 

k,< 



c" = (cf, ...,<£_!,()) (A; = 1,2). Let <r*- - ^ ,.,u n 



CTi ,.,OV, 



be defined as 



aY = J^°r^ +1 e u i' e( - r Ur (i = l,..,n). We claim that 
(3.11) a k ' € = ((Ti, .,0 + o(l), A; = 1,2. 

and 



(3.12) 



9a 



do* 
5co 



+ o(l), 



l,.,n, j = l,..,n-l, A; = 1,2. 



Assuming (I3.1ip and (13. 12f> for the moment. Now the matrix 



/ d ci a 1 



k.c 



c„-i°l 



k.c 



a 



n-l 



is non-singular at c k (k = 1, 2) for e small. On the other hand, a l,e and cx 2 ' e 
both satisfy 



(3.13) < 



f Af(ff*.«) := £ i£/ 2(2 + Pi)a k,e - IWay + dtf)*fV' e 







A} > 0, £ J £ J. 
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We use II e to represent the hyper-surface described as above. For a 2 ' 6 = 
(<T 2,e , .., an e ) G n e , we can find c 1,e = (c\ ,e , .., c^Li) such that 

cj' e = c} + o(l), j = l,2,..,n-l 

and a solution u l,t of (|3,10p with the initial condition (c\ ,e , ..,c^ l 'l 1 ,0) such 
that 

f'°° -i 

/ r ^ +1 e a ^dr = a 2 ' e , j = 1, 2, ,,n - 1. 

After using Aj(cr 2,<E ) = in (|3.13j) we have 

00 -i 

Then the difference between c 1 and c 2 implies c 1,e ^ c 2 for e small. A 
contradiction to the uniqueness property satisfied by the system (|3.10p . 

To finish the proof we now verify (|3.1ip and (|3. 12j) . Here we require 
e S (0, So) where So is so small that the matrix (oy + eSij) nxn is non-singular 
for all e G (0,S ). 

For u k , there exists Rq large such that for r > Rq and some 5 > 0, 

K fc )'(r)r < -2 - ft - 2S, i = l,...,n, k = 1,2. 

For 5o small we have u^' £ converges uniformly to over < r < Rq. For 
r = Ro we have 

(«J' e (r))V < -(2 + Pj + 8) at r = i? , < e < S . 

Then by the super-harmonicity of u k ^ e it is easy to show 

(^• e (r))V<-(2 + /3, +<5) for r > R . 
Thus, 3C > and Ri > Ro such that 
(3.14) r ft e u i' e W < c- r -(2+5) for r > Rx 

Hence for k = 1,2, 

k r°° 
eV'V^dr = / eM'V' +1 dr+o(l) = 0}+o(l), j = 1, ..,n. 



oc 







([3~TTj) is verified. To show ff3TT2|) 

(3.15) 7^-=/ r ft+1 e V (r) — i = l,..,n, fc = l,2. 
satisfies the following linearized equation: 
-A(— 3-) = V(ay + eSi^e^ i = l,..,n, Z = l,...,n-1. 
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By Lemma 12.21 

<-> k,e 

(3.16) I— i-(r)| < Clnr, r > 2, i = l,..,n, / = l,..,n-l 

dq 

where the constant C is independent of e G (0, <5o). Moreover, for any fixed 

R > 0, (r) converges uniformly to gj-(r) over < r < i? with respect 
to e. Using the decay estimates (|3.14p and (|3. 16[) in f|3. 15j) we obtain (|3. 12[) 
by elementary analysis. Proposition 13. II is proved in all cases. □ 

3.3. The proof of the second statement ofTheorem ll.il Our proof is 
based on the uniqueness result and is completely different from the method 
employed in [22]. We divide the proof into two cases according to the diag- 
onal entries of A. 
Case one: an > for all i G I. 

In this case, by Lemma 12.41 for any c±, c n _i G K, there exists a unique 
finite energy solution u = (ui, -u n ) such that Uj(0) = q for % = l,..,n — 1 
and u n (0) = 0. By Proposition 13.11 there is a bijection between the initial 
condition (ci, .., c n _i, 0) and LT2 (see the notation in the proof of Proposition 
13.11) . Thus Theorem II .11 is proved in this case. 

Case two: There exists io € I such that aj 0) j = 0. 
Let a € II2, then for e > we consider 

' -Auj = E i6 /(oii + e^lxfieW, M. 2 , 

< ul(0) = cl, i = l,...,n-l, <(0)=0, 

k f^r^ +1 e u ^dr = at, iel 
where a e = (erf, ...,<t^) is a point on the hyper-surface 

:= {a = (ai,...,a n ); o t > 0, V« G /, Af(<r) = 0, A^(o-) > 0, V0 C J C 1} 
such that a e — > a as e — >• 0. Here we recall that Aj(cr) is defined as 

The vector (cf , c^_ 1; 0) is the initial condition corresponding to cr e . 
Now we claim that 

(3.17) maxc; < C i = l,...,n — 1 

for some C > independent of e. Indeed, if this is not the case, without 
loss of generality we assume c\ is the largest among c\ and tends to infinity. 
Re-scale u e according to c\ to make the maximum of all components at 
equal to 0. The re-scaled system has to converge in Cf oc (M. 2 ) norm to 
a partial system. Indeed, the first component converges because all the 
components are bounded. The n — th component tends to —00 because 
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the initial condition is before the scaling and all components are non- 
increasing. Therefore for the limit function v = (vi,...,v n ) without loss of 
generality we assume v m +i = ...v n = for some 1 < m < n. For i = 1, m 
we easily observe that 

POO 

(3.18) Oi := / r^ i+1 e Vi{r) dr <a u i = l,...,m. 

Jo 

The reason is for each fixed R > we have 

r R 

/ r^ +1 e v ^dr < a\ + o(l), i = l,...,m. 
J o 

Clearly (vi,...,v m ) satisfies 

A «i + E7=i ayr^e^ =0, i = 1, m, 

J °° r^ i+1 e Vi ^dr < <Tj, i = 1, .., m. 
By Lemma 12.11 

m 

(3.19) s ^aij(jj>2 + (3i, i = l,...m. 

3=1 

We claim that a = (<ri, cr n ) with a m+ i = ... = a n = satisfies A/ (<r) = 0. 
Indeed, let v m+ i = ... = v n = and Hj = 1 if i = 1, ...,m and Hi = for 
i = m + 1, n. Then the system for v can be written as 

n 

Avi + ^2 ''' r ' =0, i = 1, n. 
i=i 

Apply the standard method to obtain the Pohozaev identity to the system 
above we have Aj(ct) = 0. Let J = {l,...m} we have Aj(a) = 0. Let 
Zi = &i — <Ji- From the definition of <7j we know that z% > for i = 1 , . . . , m 
Since 1 < m < n we have Aj(a) > 0, Aj(cr) — A j(a) > gives 

(3.20) Yl (E " ( 2 + ft))* + E oy^i " (2 + A)KJ < 0. 

Since X)jg j a «i^i > 2 + A for all i £ J, we also have ^ g j a?j<7j > 2 + /3j 
for all i € J because crj > crj. Clearly (|3.20p is impossible. (|3.17p is proved. 
Similarly there is a lower bound for cf, ...,c^_ 1 . As e — > 0, the ii e converges 
to u that corresponds to u. Theorem 11.11 is proved in both cases. □ 

4. Proof of Theorem 11.21 

The proof of Proposition 4.1 in [22] can be readily applied to prove The- 
orem [L2l We include it for the convenience of readers. 
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>2. 



For A > 0, let u x (xi,x 2 ) = u;(2A - x\,x 2 ). Set S A = {x G R 2 ; x\ > A.} 

1 1 ■■■■> u n) 



and T\ be the boundary of T,\. The equation for u x = (u x , ■ ■■,u x ) is 



(4.1) Au x + ^a ii |x A |e n ^ =0, iel 

jei 

where x x = (2A — x\,x 2 ). Set w x = u x — Ui to be defined in T,\ for A > 0. 
For w x we have 

Aw i +^2aij\xf l e^Wj = -y2aij(\x x \fo - \x\^)e u ^ 

3 3 

where 



Since /3, < for all i G /, 

(4.2) A^ A + ^2 aij \xf l e^w x < 0. 

j 

Let / = loglog(jx| + 3), then 

Af(x) :>) ' 



r(r + 3) 2 log(r + 3) (r + 3) 2 log 2 (r + 3) ' 
Therefore for any e > 0, there exists C(e) > such that 

(4.3) T""^ 7 ' r>C{€) - 

Let z x = w x / f , then the following lemma holds. 



Lemma 4.1. There exists R > independent of A suc/i that for A > 0, if 
xq is a point where a negative minimum of mm{z x , ..z x } is attained, then 
x G B R . 

Proof of Lemma I4.lt From (|4.2p we obtain 

(4.4) Az x + 2Vz x ^- + z i^f + <Hje$z$ < 0. 

J J 3 

Suppose z x (xo) = minj z x (xo) < and xq is where the negative minimum 
for z x is attained. Here we note that the global minimum of z x should be 
attained. Indeed, by Lemma [27TT 

Ui(x) = -mi log \x\ + Ci + 0(\x\~ 6 ) 

when \x\ is large. Thus, for A > 0, since \x x \ < \x\, 

w x (x) = u x {x) - Ui(x) > 0{\x\~ 5 ), \x\ » 1. 
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Thus ]im| a .|_ >00 zf{x) > 0. Let J = {j G /; zj(xo) < 0}. Here we observe 
that the image of the origin is not in J, because of the decay rate of U{. We 
rewrite (|4.4p as 

(4.5) Az? + 2Vz^ + + £ ayetf^ < 0. 

in a small neighborhood of xq. Then at xq, 

Az}(x )>0, 2V^ A (xo)^^ = 0. 

For j G J, since Wj(xq) < 0, we have Uj (xq) < uj(xo), so if \xq\ is large, by 
Lemma 12. 11 ^ ~ |x| _2_<5 for x close to xo and some 5 > 0. Thus 

J>,e^ ^(x ) < :,i.r„)^a ;/ ^ : ''- < Cz?(x )\x \- 2 - 5 . 
je.J j 

On the other hand if |xo| is large 

zHxo)j^>\zH*o)\\x \- 2 -*. 

Therefore by choosing e < 5/2 we see that (|4.5p can not hold if |xo| is large. 
Lemma 14.11 is established. □ 

By Lemma 14.11 and Lemma 12.11 min{u; A , w^} > in for A suffi- 
ciently large. Thus set 

A := inf{A > 0; min-Ju^, w*} > in S A }. 

Lemma 4.2. A = 0. 

Proof of Lemma 14. 2t If A > 0, we first prove that > in for 
all i G I. Indeed, let Io = {i £ I; = 0}. If Iq is not empty, the 

irreducibility of A implies all = in S^. However, not all are 0, so for 
some i G I, we have 

Awf + ^2a i:j \xfie$wj = -^c%-(|^|^ - \x\%)e"i < 0. 
jel jel 

A contradiction. 

Next we derive a contradiction to the definition of A. Let A^ tend to A from 
the left. Thus Xk > for all large k. We can assume that minjg/ w i k < 
in T,\ k because otherwise, the strong maximum principle implies w^ h > in 
T,\ k , a contradiction to the definition of A. Therefore, let x^ be where the 
minimum of minj g /w A be attained and there is ik G I such that w^ k (xk) = 
minj e / iXe s A w^ k < 0. By Lemma I3~TI Xk G Br for some R > and all k. 
Along a subsequence {xk} converges to x G such that for some io £ ^> 
•u;^ (xq) = 0. Since we have proved that > for all i G I in S\, xq G T\. 
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However, Vwi k (xk) = leads to Vw^ (xo) = 0, a contradiction to the Hopf 
Lemma. Lemma 14.21 is established. □ 

Thus we have proved A = 0, which leads to 

Ui(-xi,x 2 ) > Ui(xi,x 2 ), Van > 0, iel. 

Moving the plane from all possible directions we obtain the symmetry of U{. 
Theorem 11.21 is established. □ 

5. Uniqueness theorem on the linearized system 

In this section we prove Theorem 11.31 The following lemma describes the 
projection of u on svakO and cos k6. 

Lemma 5.1. Let (pi ) k{ r ) satisfy 

(5.1) <% k + + Y, a l3 r^e u ^Uj,k - ^i,k = 0, < r < oo 
and 

(5.2) \(/) tt k{r)\ <Cr k {l + r)- 2k , Vr > 0, k > 1. 
// there exists f = .., f n ) such that 

(5.3) //' + -ft --f i + Y <>,,>■ < 0, r > 

r r i=i 

and 

(5.4) /j(r) > 0, Vr > 0, lim fi(r)/r h = oo, lim fi(r)r k = oo. 

r— >0 r— >oo 

T/ien = 0. 

Proof of Lemma I5.lt We only need to show <^ < 0. Suppose this is not 

the case. Then because of the assumptions on the decay rates, without loss 
of generality we assume 

(5.5) Wl (r ) = = max > 0. 

/l(ro) hr fi(r) 

Note that the maximum can be attained because of the decay assumptions 
on (pi k an d (j5.4|) . The equation for w\ , after simple derivation, is 



fl r fi r fx r 2 



Near tq, wi(r) > 0. Thus in the neighborhood of ro, using (|5.3[) we have 
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The left hand side of the above is non-positive when evaluated at ro, while 
the right hand side is non-negative. A contradiction. Lemma 15.11 is estab- 
lished. □ 

Proof of Theorem II. 3t 

Let fi = —u'Ar). Direct computation shows that 

3 j 

Since all ft < 0. / = {fx, ...,f n ) satisfies ([OD and (JED for all k > 2. Let 
(fi k = ...,</>*) be the radial part of the projection onto, say, sinkO. Then 
(ft k satisfies (|5.ip . Since fc is bounded, it is easy to apply standard ODE 
theorem to obtain that (|5.2p also holds. Thus all the projections on sin k9 
and cos kO are all zero for f3 < 0. 

Finally we prove that for the projection on sin# or cos 6 is also zero. Let 
1 = (01, i) -j 01, n) be the projection of (ft on sin#. Then we have 

<f>'U + Vm - ^<t>U ■ Y."W !, '" J OLj = 0. 

j 

Since (ft 1 is bounded, the standard ODE theory implies that (ftn behaves like 
0(l/r) at infinity and like 0(r) near 0. We shall use / = (— u[, —u' n ) as 
the function to majorize (ft 1 . To apply the same argument as in the proof 
of Lemma 15. 1\ The problem is the maximum may tend to or infinity. We 
first prove that this can not happen at 0: 

(5.6) <fti/ fi is strictly increasing near if (fti is positive near 0. 

Clearly once (]5.6p is proved, (ft 1 = 0, thus Theorem 11.31 would be established. 

Now we prove (15. 6ft . Let zi = (fti^/r and Fi = fi/r. Direct computation 
yields 

z" + -z[ + aijr^e Uj zj = 0, r > 0. 
j 

and 

i j 

Since is positive near 0, Zj(0) > (if Zi(0) = 0, there is no need to 
consider this case, as the maximum can not tend to 0). Easy to see that 
near 0, 

Zi (r) = z t (0)+Y,O(r^ +2 ) 

j 

and 
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Proving 4>i t i/fi to be increasing near is the same as proving that Zi/Fi is 
increasing near 0. Since < 0, one immediately sees that Zi/Fi is increasing 
near 0. 

Next we prove that Zi/Fi is decreasing if Z{ is positive at infinity. Assume 
z i( r ) = Qi/f 2 + 0(r 3 ) at infinity. We have known that, for some Si > 0, 

u i{ r ) = ~ m i l°g r + °i + 0(r~ Si ) r > 1. 

Thus 

e n 8 (r) = e c lr - mi + (r- m <~ 5 >), r » 1. 

We obtain, by integration on the equation for Zi, that 



.7 



Correspondingly to compute F{, we use the equation for m to obtain 

(ruj(r))' = -J2 a ij rf>j+leUi = -J2 a v rl+ ^~ mie ° j +0(r 1+ h- m i- s i). 
j j 

Using ru^(r) — > —rrii at infinity, we have 



g j 

ru'i(r) = -rrn + V a„ _ r 2 +ft-^ + 0(r 

m j Pj * 

Consequently 



2+/3j—mj—Sj •• 



— = 5 -E- g4^ to(^). 



r r 2 — j3j — 2 



j- Pj ^ 

Our goal is 

— Zji 7 / < near infinity 

when Zi is positive near infinity. Using the expressions above it is enough to 
show if the following is negative: 

(5-7) E^Ha + ^V*)^" mj " 2 - 

j "v Pj 

By qi/rrii = maxjgj qj/mj, qi > and < for all i £ I, we have (|5.7p . 
Therefore Zj/-Fj is decreasing near infinity. Theorem 11.31 is proved. □ 
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6. Appendix 

In this appendix, we prove the ODE lemmas stated in section two. 
Proof of Lemma I2.lt The proof is standard ( for example, see |20j). We 
include it for the convenience of the reader. Let 

(6.1) Wi (x)= / (-—] og \x-y\ + — ]og(l + \y\))^ai j \yf'e^Wdy. 
J R 2 2tt 2tt ^— ' 

Clearly wi is well defined and satisfies 

-Aw i (x) = Y J a ij\x\' 3j e u > {x \ M 2 
j 

and 

A( Ui -Wi) = 0, M 2 . 
By Lemma 4.1 in |40j m < C onl 2 . Next we claim that 

,. Wi{x) 

lim j — - = —mi. 

|a;|-K» log \x\ 

To see the above, it is easy to obtain for e > 0, there exist R(e) >> 1 and 
R\ >> R such that for Ixl > R\ 



1 f -log|s-y| + log(l + \y\) v-^ I \Bj Ui |^ 
2vr J B \og\x\ l3m 11 " 



Also 



Jm?\b r log|x| ^ 

Thus Ui — Wi < Clog(l + \x\), which leads to 

(6.2) Ui = Wi + Ci 
for some Cj £ R. Next we claim that 

(6.3) mi - Pi > 2 for all i € I. 
Indeed, if this is not the case, there exists io G I such that 

m io ~ A = min{m! - f} u ..,m n - /3 n } < 2. 
By (J6ZEJ) and ([El we have 

Ui (x) = ^- [ {-\og\x-y\+\og(l + \y\))Y^a i0J \y\^e u ^dy-C. 
Easy to check 

- log \x-y\+ log(l + \y\) > - log(l + |x|), 

thus 

u io (x) > -mi,, log(|x| + 1) - C > -(2 + /3 io )log(|x| + l)-C 
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a contradiction to J R2 {x^o e Ui o( x > < oo. (|6.3p is established. Now Ui(x) can 
be written as 

(6.4) Ui (x) = -— I \og\x-y\ y^Oijlyf^e^dy + Ci. 

Q can be determined as in 
0(r _<5 ). To see this we set 



3 

Ci can be determined as in the statement. Finally we derive the error term 
-6\ 



Ei = {y, \y\<\x\/2.} E 2 = {y; \ y - x \<^.} E 3 = R 2 \ (E 1 U E 2 ). 
Using e Ui ^ = 0(\y\~ mi ) in (I6.4p one obtains 

-— / log\x-y\y^a ij \y\^e u ^ y) dy = -m i \og\x\ + 0(\x\- s ). 

Similarly by elementary estimates 

/ log \x - y\Y^ aij \y\^e u ^dy = 0(\x\- 5 ). 
Je 2 ue 3 ■ 

The gradient estimate for u\ is obtained by standard estimates. Lemma 12.11 
is established. □ 

Proof of Lemma 12. 2\ Let ip(t) = {ipi{t), ...,i/j n (t)) be defined as 

ipi(t) = &(e*), i E /. 

Then if) satisfies 

<(*) + aije Viie * )+{7+Pi Hj(t) = 0, -oo < t < oo, is I. 
j 

Let ip n+ i = ip'i,..., i>2n = tp'n and F = (ipi, .., V>2n) T , then F satisfies 

F' = MF 

where M = ^ ^ jj.Bisanxn matrix with By = —aije Uj( - et ^ +( - 2+l3 ^ t . 
For t > 1, the solution for F is 

(6.5) F(t) = lim e eM ^...e eM ^F(0). 

N— >oo 

where to, •••,tjv satisfy tj = je, j = 0, .., N, e = t/N. Since Uj(e*) + (2+/3j)t ~ 
(— mj + 2 + /3j)t when t is large and > 2 + ft, we have ||B|| ~ e _<5t for 
some 5 > and t large. With this property we further have 

(6.6) ||M|| fc < Ce~ k5lt , k = 2, 3, ... t > 
for some 5i > 0. Using (16. 6j) in (16.5 we have 

||F(t) || = 0(t), t>l. 
Lemma 12.21 is established □ 



(fc+1), 
u) ( r 
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Proof of Lemma I2.3t If a solution u = (u±, ...,u n ) exists, it would satisfy 

u,(r) = c i (0)-Va !J t^ +1 (logr-logt)e u ^dt, i = l,...,n. 
3=1 Jo 

We first prove the existence of a solution on < r < 5 for some small 5 > 
by iteration: Let = (0, ..,0) and 

) = Ci(0)- Voy /V' +1 (logr-logi) e u fWcft, i = l,...,n. 
i=i 7o 

For (5 > small and r G (0, (5), since + 1 > —1, it is easy to see that such 
a sequence converges. Therefore the existence of a solution for over (0, 5) 
is proved. The existence for r G (<5, oo) clearly holds because of the right 
hand side is a Lipschitz function of it. The proof of the uniqueness of the 
solution is the same as that in Lemma 12.51 later in the section. Lemma 12.31 
is established. □ 

Proof of Lemma 12. 4t By Lemma 12.31 a solution to (|2.2p exists for r > 0. 
We just need to show that 

POD 

/ e Ul(r) r ft+1 dr < oo, Vi G /. 

Let Vi(t) = Ui(e ) + (2 + j3i)t (i G /), then v = {v\, ...,v n ) satisfies 
v" (t) + ^ aije V]{i) = 0, -oo < i < oo, i £ I. 

3 

From the equation for u we have 

ru' l (r) = - [ Vo !3 e^ (s »sft +1 (is<0, r > 0, i G I. 
Jo j 

The last inequality is strict because aij > and not all equal to 0. Conse- 
quently v'^t) < 2 + ^ for t G E. Fix i G R we have, for f > i > 

^(t) = ^(t„) _ /" J^Oij-e^Wds < v'So) - an f e v ^ s) ds, i G /. 

J to j J to 

Since an > there exists t > to such that v^it) < 0. Choose t^(ii) = — 5 < 
for some 5 > 0, then we see 

which is equivalent to u,(r) < (—2 — — 5) log r + C for r > e* 1 . Therefore 
J °° e u ii r ) r^ i+1 dr < oo. Lemma I2TH is established. □ 

Proof of Lemma 12. 5t The proof is standard, we include it for the conve- 
nience of the reader. Clearly we only need to show that fa = in (0, 6) for 
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5 > small. Write <fii(r) as 

Mr) = ~ f~ f (yW ft+ W%(i)di)d S 

rr 

= - Va^+W^i) (log r- log i)tft 
Jo jei 

Let a = min{/?i, f3 n } + 1, since all > —2 we have a — e > — 1 for some 
e > small. For the e we choose <5 > small so that log r — log t < t~ e for 
r < 5 and t < r. Thus 

|^i(r)|<c/'V- e y;|^(t)|dt J r<5 
Jo j 

for some C. Let 0(r) = Eie/ \M r )\ and = /J" ^I^C*)!*' then 

F'(r) - Cr a - e F(r) < 0, F > 0, F(0) = 0. 
Since a — e > —1, i 7 = 0. Lemma 12.51 is established. □ 

References 

[1] D. Bartolucci, Compactness result for periodic multivortices in the electroweak theory, 

Nonlinear Anal., 53 (2003), 277297. 
[2] D. Bartolucci and G. Tarantello, The Liouville equation with singular data: A 

concentration- compactness principle via a local representation formula, J. Differ- 
ential Equations, 185 (2002), 161180. 
[3] D. Bartolucci and G. Tarantello, Liouville type equations with singular data and 

their appli- cations to periodic multivortices for the electroweak theory, Comm. Math. 

Phys., 229 (2002), 347. 
[4] D. Bartolucci, C. C. Chen, C. S. Lin and G. Tarantello, Profile of blow-up solutions 

to mean field equations with singular data, Comm. Partial Differential Equations, 29 

(2004), 1241 1265. 

[5] W. H. Bennet, Magnetically self-focusing streams, Phys. Rev. 45 (1934), 890897. 
[6] H. Brezis and F. Merle, Uniform estimates and blow-up behavior for solutions of 

-Au = V(x)e u in two dimensions, Comm. Partial Differential Equation, 16 (1991), 

12231253. 

[7] L. Caffarelli and Y. Yang, Vortex condensation in the Chern-Simons Higgs model: an 

exis- tence theorem, Comm. Math. Phys., 168 (1995), 321336. 
[8] E. Caglioti, P. L. Lions, C. Marchioro and M. Pulvirenti, A special class of station- 
ary flows for two-dimensional Euler equations: A statistical mechanics description, 

Comm. Math. Phys., 143 (1992), 501525. 
[9] E. Caglioti, P. L. Lions, C. Marchioro and M. Pulvirenti, A special class of stationary 

flows for two-dimensional Euler equations: A statistical mechanics description, part 

II, Comm. Math. Phys., 174 (1995), 229260. 
[10] S. Chanillo and M. Kiessling, Rotational symmetry of solutions of some nonlinear 

problems in statistical mechanics and in geometry, Comm. Math. Phys., 160 (1994), 

217238. 

[11] S. Chanillo, M. K-H Kiessling, Conformally invariant systems of nonlinear PDE of 

Liouville type. Geom. Funct. Anal. 5 (1995), no. 6, 924-947. 
[12] S. Y. Chang, M. Gursky and P. Yang, The scalar curvature equation on 2- and 3- 

spheres, Calc. Var. and PDE, 1 (1993), 205229. 



24 CHANG-SHOU LIN AND LEI ZHANG 

[13] S. Y. Chang and P. Yang, Prescribing Gaussian curvatuare on S2, Acta Math., 159 
(1987), 215259. 

[14] C. C. Chen and C. S. Lin, Estimate of the conformal scalar curvature equation via 

the method of moving planes. II, J. Differential Geom., 49 (1998), 115178. 
[15] C. C. Chen, C. S. Lin, Sharp estimates for solutions of multi-bubbles in compact 

Riemann surfaces, Comm. Pure Appl. Math. 55 (2002), no. 6, 728-771. 
[16] C. C. Chen, C. S. Lin, Topological degree for a mean field equation on Riemann 

surfaces, Comm. Pure Appl. Math. 56 (2003), no. 12, 1667-1727. 
[17] C. C. Chen, C. S. Lin, Mean field equations of Liouville type with singular data: 

Sharper estimates, Discrete and continuous dynamic systems, Vol 28, No 3, 2010, 

1237-1272. 

[18] C. C. Chen, C. S. Lin and G. Wang, Concentration phenomena of two- vortex solutions 

in a Chern-Simons model, Ann. Sc. Norm. Super. Pisa CI. Sci. (5), 3 (2004), 367397. 
[19] W. Chen and C. Li, Classification of solutions of some nonlinear elliptic equations, 

Duke Math. J., 63 (1991), 615622. 
[20] W. X. Chen, C. M. Li, Qualitative properties of solutions to some nonlinear elliptic 

equations in R 2 . Duke Math. J. 71 (1993), no. 2, 427439. 
[21] S. Childress and J. K. Percus, Nonlinear aspects of Chemotaxis, Math. Biosci. 56 

(1981), 217237. 

[22] M. Chipot, I. Shafrir, G. Wolansky, On the solutions of Liouville systems. J. Differ- 
ential Equations 140 (1997), no. 1, 59-105. 

[23] M. Chipot, I. Shafrir, G. Wolansky, Erratum: "On the solutions of Liouville sys- 
tems" [J. Differential Equations 140 (1997), no. 1, 59-105; MR1473855 (98j:35053)]. 
J. Differential Equations 178 (2002), no. 2, 630. 

[24] P. Debye and E. Huckel, Zur Theorie der Electrolyte, Phys. Zft 24 (1923), 305325. 

[25] G. Dunne, Self-dual Chern-Simons Theories, Lecture Notes in Physics, vol. m36, 
Berlin: Springer- Verlag, 1995. 

[26] J. Hong, Y. Kim and P. Y. Pac, Multivortex solutions of the Abelian Chern-Simons- 
Higgs theory, Phys. Rev. Letter, 64 (1990), 22302233. 

[27] R. Jackiw and E. J. Weinberg, Selfdual Chern Simons vortices, Phys. Rev. Lett., 64 
(1990), 22342237. 

[28] J. Jost, C. S. Lin and G. Wang, Analytic aspects of the Toda system. II. Bubbling 
behavior and existence of solutions, Comm. Pure Appl. Math., 59 (2006), 526558. 

[29] J. Jost, G. Wang, Classification of solutions of a Toda system in R 2 . Int. Math. Res. 
Not. 2002, no. 6, 277-290. 

[30] J. Jost, G. Wang, Analytic aspects of the Toda system. I. A Moser-Trudinger inequal- 
ity. Comm. Pure Appl. Math. 54 (2001), no. 11, 1289-1319. 

[31] J. Kazdan and F. Warner, Curvature functions for compact 2-manifolds, Ann. of 
Math., 99 (1974), 1447. 

[32] E. F. Keller and L. A. Segel, Traveling bands of Chemotactic Bacteria: A theoretical 
analysis, J. Theor. Biol. 30 (1971), 235248. 

[33] M. K.-H. Kiessling and J. L. Lebowitz, Dissipative stationary Plasmas: Kinetic Mod- 
eling Bennet Pinch, and generalizations, Phys. Plasmas 1 (1994), 18411849. 

[34] Y. Y. Li, Harnack type inequality: the method of moving planes, Comm. Math. Phys. 
200 (1999), no. 2, 421-444. 

[35] C. S. Lin, An expository survey on the recent development of mean field equations, 
Discrete Contin. Dyn. Syst., 19 (2007), 387410. 

[36] C. S. Lin and M. Lucia, Uniqueness of solutions for a mean field equation on torus, 
J. Differential Equations, 229 (2006), 172185. 

[37] C. S. Lin and M. Lucia, One-dimensional symmetry of periodic minimizcrs for a mean 
field equation, Ann. Sc. Norm. Super. Pisa CI. Sci. (5), 6 (2007), 269290. 

[38] C. S. Lin, J. C. Wei, D. Ye, Classifcation and nondegcncracy of SU(n + 1) Toda 
system, preprint. 



LIOUVILLE SYSTEM WITH SINGULARITY 



25 



[39] C. S. Lin, J. C. Wei, C. Zhao, Sharp estimates for fully bubbling solutions of a SU(3) 

Toda system. Preprint. 
[40] C. S. Lin, L. Zhang, Profile of bubbling solutions to a Liouville system. (English, 

French summary) Ann. Inst. H. Poincar Anal. Non Linaire 27 (2010), no. 1, 117-143. 
[41] C. S. Lin, L. Zhang, A topological degree counting for some Liouville systems of mean 

field equations, Comm. Pure Appl. Math, volume 64, Issue 4, pages 556-590, April 

2011. 

[42] M. S. Mock, Asymptotic behavior of solutions of transport equations for semiconduc- 
tor devices, J. Math. Anal. Appl. 49 (1975), 215225. 
[43] M. Nolasco and G. Tarantello, On a sharp type inequality on two dimensional compact 

man- ifolds, Arch. Rational Mech. Anal., 145 (1998), 161195. 
[44] M. Nolasco and G. Tarantello, Double vortex condensates in the Chern-Simons-Higgs 

theory, Calc. Var. and PDE, 9 (1999), 3194. 
[45] M. Nolasco and G. Tarantello, Vortex condensates for the SU(3) Chern-Simons theory, 

Comm. Math. Phys., 213 (2000), 599639. 
[46] J. Prajapat and G. Tarantello, On a class of elliptic problems in R2: Symmetry and 

unique- ness results, Proc. Roy. Soc. Edinburgh Sect. A, 131 (2001), 967985. 
[47] T. Ricciardi and G. Tarantello, Vortices in the Maxwell-Chern-Simons theory, Comm. 

Pure Appl. Math., 53 (2000), 811851. 
[48] I. Rubinstein, Electro diffusion of Ions, SIAM, Stud. Appl. Math. 11 (1990). 
[49] J. Spruck and Y. Yang, Topological solutions in the self-dual Chern-Simons theory: 

Existence and approximation, Ann. Inst. H. Poincare Anal. Non Lineaire, 12 (1995), 

7597. 

[50] J. Spruck and Y. Yang, On multivortices in the electroweak theory I: Existence of 

periodic solutions, Comm. Math. Phys., 144 (1992), 116. 
[51] J. Spruck and Y. Yang, On multivortices in the electroweak theory II: Existence of 

Bogo- molnyi solutions in R2, Comm. Math. Phys., 144 (1992), 215234. 
[52] G. Tarantello, A Harnack inequality for Liouville-type equations with singular 

sources, Indiana Univ. Math. J., 54 (2005), 599615. 
[53] G. Tarantello, A quantization property for blow up solutions of singular Liouville-type 

equa- tions, J. Funct. Anal., 219 (2005), 368399. 
[54] G. Tarantello, Analytical aspects of Liouville-type equations with singular sources, 

Stationary partial differential equations, Vol. I, 491592, Handb. Differ. Equ., North- 
Holland, Amsterdam, 2004. 
[55] G. Tarantello, Selfdual gauge field vortice. An analytical approach, in Contributions 

in Nonlinear Differential Equations and their Applications, 72, Birkhauser, Boston, 

2008. 

[56] Y. Yang, Solitons in Field Theory and Nonlinear Analysis, Springer- Verlag, New 
York, 2001. 

[57] L. Zhang, Blow up solutions of some nonlinear elliptic equations involving exponential 
non- linearities, Comm. Math. Phys., 268 (2006), 105133. 

[58] L. Zhang, Asymptotic behavior of blowup solutions for elliptic equations with expo- 
nential nonlinearity and singular data, Communications in Contemporary mathemat- 
ics, vol 11, issue 3, June 2009, 395-411. 

Taida Institute of Mathematical Sciences, and Center for Advanced Study 
in Theoretical Sciences, National Taiwan University, Taipei 106, Taiwan 
E-mail address: cslin@math.ntu.edu.tw 

Department of Mathematics, University of Florida, 358 Little Hall P.O.Box 
118105, Gainesville FL 32611-8105 
E-mail address: leizhang@ufl.edu 



